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Abstract. Collapsing collisionless particle systems form gravitational bound halos 
with cuspy density profiles. Also hierarchical merging of these systems produce 

04 ' remnants with cuspy central density profiles. These results lead to the assumption of 

, cuspy NFW [5] profiles for the density distribution in dark matter halos. However, 

■ observed rotation curves in disk galaxies suggest dark matter halos with isothermal 
core. The same kind of problem can be found for globular clusters, which show 

^ |. cores in their density profiles but should be cuspy if they where formed through a 

' cold collapse. [4, 3] We are showing how small, massive, and compact objects can 

I efficiently transform cuspy stellar systems into density profiles with an isothermal 

o ' 
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^ ■ 2.1. Initial conditions 

O I For the realizations of cuspy systems we are choosing Dehnen models 

' with isotropic velocity dispersion. The Dehnen family of models are 

. characterized by a parameter 7 that measures the degree of central 



concentration. [2] The density profile is 



, , (3-7)M a 
C3 I 47r rT (r + a) ^ 

M being the total mass and a being the scale length. Both M and 
a are set to unity. We chose centrally condensed Dehnen models with 
inner slope 7 = 1, 3/2, and 2, yielding a density profile similar to those 
of elliptical galaxies with dense stellar nuclei. The black holes and the 
stars are initially distributed according to the same profile. 

We are creating several realizations of these Dehnen models with a 
total particle number of 9910. These realizations contain 10 black holes 
with a combined mass of 0.01 which is 1% of the total mass of the 
system. 
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Figure 1 . Comparison of the initial radial density distribution of the stars with the 
average radial density distribution of the stars after 600 iV-body time units. 



2.2. Method 

We are using NB0DY6++ for the orbit integration. NB0DY6++ is a 

high precision fourth order direct force integrator [6, 1]. Since it can 
regularize the equations of motion for particles having close encounters 
no softening needs to be applied. This is why we are able to follow the 
formation of black hole - black hole binaries transfering part of their 
binding energy into the stellar system. 

We stopped the simulations once the number of black holes bound 
to the stellar system dropped below three. This process usually took a 
few hundred dynamical times. 

2.3. Results 

Figures 1 and 2 summarize our results. Both plots compare the initial 
model with the situation after 600 A^-body time units. This corrosponds 
to an evolution over roughly 100 dynamical timescales. We have aver- 
aged the results of four realizations in order to find p{r) and (T(r) in 
Figures 1 and 2. 

Figure 1 shows the initial radial density distribution compared with 
the averaged radial density distribtion after 600 AT-body time units. 
While the impact of the black holes removes mass up to radius rcrit 
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Figure 2. Comparison of the initial 3D velocity dispersion of the stars with the 
average 3D velocity dispersion of the stars after 600 iV-body time units. 



Table I. The amount of relocated 
mass after 600 A'^-body time units 

7 rcrit Mrel 



3/2 



0.5 0.06 
0.4 0.07 
0.3 0.07 



of 0.5 which is half way to the initial scale radius for 7 = 1.0 and 
7 = 1.5, the effect is seemingly smaller in the models with 7 = 2.0. As 
Table I shows, the amount of mass relocated by the black holes remains 
constant. 

Figure 2 shows the initial 3D velocity dispersion in the models and 
the 3D velocity dispersion after 600 iV-body time units. Since models 
with 7 = 2.0 are isothermal, the evolution is quite modest for such 
models. In the case of smaller 7, Dehnen models have a zero velocity 
dispersion in the center. Our models with smaller 7 show a strong 
evolution into isothermal systems through the influence of the black 
holes. 
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3. Discussion 

We have shown that 10 black holes are capable of efficiently destroying 
a stellar cusp within a few ten dynamical timescales. In our simulations 
each black hole removes an equivalent of seven times its own mass from 
the central regions of the sytem. The resulting systems show isothermal 
cores. 
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